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5 i ABSTRACT 

q ■ We present 0.8— 5.2/im spectroscopy of the compact source at the base of a 

variable nebula (McNeil's Nebula Object) in the Lynds 1630 dark cloud that went 
into outburst in late 2003. The spectrum of this object reveals an extremely red 
continuum, CO bands at 2.3—2.5 /im in emission, a deep 3.0 /im ice absorption 
feature, and a solid state CO absorption feature at 4.7 /mi. In addition, emission 
lines of H, Ca II, Mg I, and Na I are present. The Paschen lines exhibit P Cygni 
profiles, as do two lines of He I, although the emission features are very weak in 
the latter. The Brackett lines, however, are seen to be purely in emission. The 
P Cygni profiles clearly indicate that mass outflow is occurring in a wind with a 
velocity of ~ 400 km s _1 . The H line ratios do not yield consistent estimates of the 
reddening, nor do they agree with the extinction estimated from the ice feature 
(Ay ~ 11). We propose that these lines are optically thick and are produced in 
a dense, ionized wind. The near-infrared spectrum does not appear similar to 
any known FUor or EXor object. However, all evidence suggests that McNeil's 
Nebula Object is a heavily-embedded low-mass Class I protostar, surrounded by 
a disk, whose brightening is due to a recent accretion event. 
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Subject headings: stars: variables: other - stars: pre-main sequence - stars: 
formation 



1. Introduction 

The object 2MASS J05461313-0006048 is located in the Lynds 1630 dark cloud of the 
Orion B molecular cloud/star formation complex. This very red object is associated with 
the IRAS source IRAS 05436—0007, the millimeter continuum source LMZ 12 (Lis, Menten, 
& Zylka 1999), and the submillimeter continuum source OriB55smm (Mitchell et al. 2001; 
Johnstone et al. 2001). This object appears point-like in the deep /—band images of the 
region obtained by Eisloffel & Mundt (1997). Johnstone et al. (2001) estimated a size of 
R cS ~ 6.8 x 10 3 AU and a mass of O.68M for the submillimeter source. Based on the 
IRAS and millimeter data, Lis et al. (1999) estimated a bolometric luminosity for LMZ 12 
of 2.7L Q and identified it as a young, heavily embedded Class source. This interpretation 
seemed to agree with the very red colors derived from the 2MASS colors, the fact that the 
source is invisible on the POSS plates, and the identification of the object by Eisloffel & 
Mundt (1997) as the source driving the HH23 outflow. However, Lis et al. (1999) also noted 
that this classification was not completely consistent with the characteristics of other Class 
sources, as LMZ 12 appeared to be very compact and lacked any sign of a molecular outflow. 

This object recently experienced a dramatic brightening, associated with the sudden 
appearance of a conical-shaped nebula in optical images of the region (McNeil 2004). As 
the 2MASS source sits at the base of this reflection nebula, we shall henceforth refer to it 
as McNeil's Nebula Object (MNO). Reipurth & Aspin (2004) have presented optical and 
near-infrared images of the source and its associated nebula, as well as optical and K-band 
spectra of MNO itself. Their spectra reveal strong 12 CO emission, Br 7 in emission, and a P 
Cygni profile for the Ha line. They suggest that the brightening of MNO may be related to a 
FUor- or EXor-like event. Briceho et al. (2004) have presented /—band images of MNO and 
its reflection nebula before and during the recent outburst. From optical spectra of MNO and 
a separate region in the nebula, obtained during the outburst, they deduce a spectral type 
of early B for MNO and a bolometric luminosity of ~ 219L . Andrews, Rothberg, & Simon 
(2004) have presented 10/im spectra, 12 CO (2—1), 450 /im, and 850 fim observations of MNO 
in outburst. Integrating under the spectral energy distributions (SEDs) generated from all 
available pre- and post-outburst measurements, Andrews et al. (2004) estimate that the 
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bolo metric luminosity of MNO increased by at least an order of magnitude, to ~ 30 — 9OL , 
during the outburst. Based on the shapes of the pre- and post-outburst SEDs, Abraham 
et al. (2004) and Andrews et al. (2004) conclude that MNO is a Class I/Class II transition 
object. 

Here we present near-infrared (0.8— 5.2/zm) medium resolution spectroscopy of MNO 
obtained during the outburst, approximately one month after the observations of Reipurth 
& Aspin (2004) and 1-2 days before the 10/im and submillimeter observations of Andrews 
et al. (2004). Based on these data, we suggest that MNO is a low-mass late-stage Class I 
protostar with mass outflow occurring in a dense, ionized wind. In this paper we assume a 
distance to MNO of 400 pc (Antony-Twarog 1982). 

2. Observations and Data Reduction 

MNO was observed at the NASA Infrared Telescope Facility on Mauna Kea on 2004 
March 09, 06h UT with SpeX, the facility near-infrared medium-resolution cross-dispersed 
spectrograph (Rayner et al. 2003). Ten individual exposures of MNO, each lasting 120 s, 
were obtained using the short wavelength cross-dispersed mode (SXD, covering 0.8—2.5 /xm) 
of SpeX. These were followed by ten exposures, each lasting 30 s, with the long wavelength 
cross-dispersed mode (LXD2.3, covering 2.3— 5.5/zm). The observations were acquired in 
"pair-mode", in which the object is observed at two separate positions along the slit. We 
used the O'.'3-wide slit, which yields nominal resolving powers of 2000 and 2500 for the SXD 
and LXD2.3 modes, respectively. The slit was set to the parallactic angle for the SXD 
observations. Observations of an A0V star, used as a "telluric standard", were obtained 
immediately preceding or following the observations of MNO. A set of internal flat fields 
and arc frames (for wavelength calibration) were also obtained for each mode. The airmass 
varied between 1.1 and 1.4 during the observations; however the airmass differences between 
the object and the standard were less than 0.06. The seeing was estimated to be ~ 1" at 2.2 
/iin and conditions were clear. 

The data were reduced using the Spextool package, which performs flat-fielding, non- 
linearity corrections, pair subtraction, aperture definition, optimal spectral extraction, and 
wavelength calibration of SpeX data semi-automatically (see Cushing et al. 2004). Flux 
calibration of the spectra was carried out using the method and software described by Vacca 
et al. (2003); the flux calibration is expected to be accurate to about 10%. After flux 
calibration, the individual SXD and LXD2.3 spectral orders were spliced together. The final 
spectrum spanning 0.8—5.2 /xm is shown in Fig. 1. An enlargement of the 1 — 1.3 and 1.9—2.5 
fj,m regions, after normalization of the continuum with a low order polynomial, is also shown. 
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The signal-to-noise ratio (S/N) varies across the spectral range, but is of the order of 200 
and 150 at the centers of the J and K bands, respectively. At the longest wavelengths the 
S/N of the MNO spectrum is limited by that of the AOV telluric standard spectrum (S/N 
~ 40). 

3. Results 

The 0.8—5.2 /iin spectrum of MNO exhibits a very red continuum with a number of 
emission and absorption features, the most prominent of which are the broad ice absorp- 
tion feature at 3.1 pm, the deep CO solid-state absorption feature at 4.7 /zm, the CO ro- 
vibrational Au = +2 transitions at 2.3— 2.4/im in emission, and the Pa a, Pa /3, Br a, and 
Br 7 emission lines. Also present are weaker emission features of Ca II at 0.85 — 0.86 /zm, 
Fe I at 1.16 - 1.20 /im, Mg I at 1.50 /im, Na I at 2.21 /im, and Ca I at 2.26 //m. Although 
the Brackett lines appear to be purely in emission, the Pa f3, 7, and 5 lines have P Cygni 
profiles, which provides clear evidence of mass outflow. Close inspection of two He I lines 
at 1.083 jum and 2.058 /im also reveals P Cygni profiles, with weak emission components. 
From the blue edges of the P Cygni profiles of the Paschen and He I lines we estimate a 
maximum outflow velocity of ~ —550 km/s; the absorption is centered at ~ —400 km/s. A 
weak emission line at 2.12 /im due to H 2 may also be present, which would provide some 
evidence for shocked gas near MNO. However, we see no other signatures of H 2 emission in 
our SpeX data. 

The estimated magnitudes in the 2MASS photometric system at an airmass of one, as 
synthesized from our spectrum, are J = 10.90, H = 8.85 and K s = 7.30, with estimated 
uncertainties of about 10%. The estimated magnitudes in the Mauna Kea Observatory 
photometric system (Simons & Tokunaga 2002; Tokunaga, Simons, & Vacca 2002) at an 
airmass of one are J = 10.83, H = 9.00, K = 7.20, V = 4.98 and M' = 4.10. These values 
are in reasonable agreement with the magnitudes estimated by Reipurth & Aspin (2004) and 
suggest that the absolute flux calibration of our SpeX data is reliable. Comparison with the 
pre-outburst 2MASS magnitudes indicates that the source brightened by 3.84, 3.31, and 2.97 
mags in J, H, and K s respectively, with an associated color change of A(J — H) = —0.53 
and A(H — K) — —0.34 mags. Hence, in the NIR, MNO was bluer during outburst than in 
quiescence. 

The depth of the 3/im ice feature implies an optical depth r(3.0) ~ 0.7. From the 
correlation between r(3.0) and Ay given by Whittet et al. (1988), we infer a reddening of 
Ay ~ 11. We have used this value and the reddening law of Rieke & Lebofsky (1985) 
to correct the observed emission line fluxes for reddening. The observed fluxes, equivalent 
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widths (EWs) and dereddened luminosities of many of the strongest features seen in our 
spectrum of MNO are given in Table 1. 

The integrated flux in the 0.8—5.2 /im spectral range corresponds to a luminosity of 
~ 9.2L Q before any reddening corrections. This is a factor of ~ 3 larger than the quiescent 
bolometric luminosity estimated by Lis et al (1999) from the combination of IRAS and 
millimeter data, and a factor of ~ 2 larger than that estimated by Abraham et al. (2004) 
and Andrews et al. (2004) from integrating under the entire pre-outburst SED. Applying an 
extinction correction of Ay = 11, we find an intrinsic luminosity over this wavelength range 
of ~ 27 L Q . 

Assuming that a ~ 5OL increase in the bolometric luminosity during the outburst 
(Andrews et al. 2004) is due to accretion onto a O.5M pre-main sequence star and adopting 
the relation between radius and mass given by Kenyon, Calvet, & Hartmann (1993), we 
estimate a mass accretion rate for MNO of ~ lO _5 M yr _1 during the outburst. Using the 
outflow velocity derived from the P Cygni profiles, and assuming that HH23A represents a 
termination point for the outflow, we calculate a dynamical timescale for MNO of ~ 800 yr, 
significantly smaller than the value estimated by Eisloffel & Mundt (1997) using tangential 
velocities derived from proper motions. 



4. Discussion 

Reipurth & Aspin (2004) and Briceno et al. (2004) have suggested that MNO is a FUor 
or EXor object. However, our NIR spectrum of MNO does not resemble any existing NIR 
spectra of either type of object. Although a K— band spectrum of EX Lupi during outburst 
exhibited weak Br a in emission, the Na I, Ca I and the CO ro-vibrational Av = +2 transi- 
tions were in absorption (Herbig et al. 2001). In addition, the Balmer lines exhibited reverse 
P Cygni profiles (Lehmann, Reipurth, & Brandner 1995). Similarly, the 12 CO transitions 
were in absorption in the K— band spectra of four known FUors (Reipurth & Aspin 1997). 
While the brightness increase in MNO is larger than that for an EXor, our estimated accre- 
tion rate for MNO is well below the typical values found for FUors (Hartmann & Kenyon 
1996). Nevertheless, there are clear indications that MNO is a young stellar object. 

When the differences in reddening are taken into account, our 1.9— 4.1/xm spectrum of 
MNO is seen to be nearly identical to that of the Class I protostar SSV59 presented by 
Simon et al. (2004), although the emission features are somewhat weaker in the spectrum of 
SSV59. This supports the suggestion, based on the low pre-outburst bolometric luminosity 
(2.7-4.2 L ; Lis et al. 1999; Abraham et al. 2004; Andrews et al. 2004) and the submillimeter 
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continuum slope (Andrews et al. 2004), that MNO is a late-stage Class I object. The fact 
that SSV59 also illuminates a reflection nebula and has a bolometric luminosity similar to 
that of MNO strengthens this conclusion. 

The origin of H lines in the spectra of Class I sources is controversial. In models 
developed by Natta et al. (1988) and Hartmann & Kenyon (1990) the H lines are generated 
in a magnetically driven, outflowing wind. However, based on the symmetric, purely emission 
profiles of Br7 in the spectra of such sources, Najita et al. (1996) argued that the H lines are 
produced by mass infall, rather than outflow. Furthermore, Muzerolle, Hartmann, & Calvet 
(1998) used the Br7 luminosity from Class I sources as a measure of the accretion luminosity 
and hence the mass accretion rate in these objects. However, neither the current magnetically 
driven wind models nor the magnetospheric accretion models are able to reproduce all of the 
various features of the H line profiles seen in the spectra of Class I sources (see Nisini et al. 
2004 for a discussion). Nisini et al. (2004) have recently argued that the Br line strengths in 
the NIR spectrum of HH 100 result from an accelerating, ionized wind. 

The Paschen and two He I lines seen in our spectrum of MNO exhibit unmistakable P 
Cygni profiles. The strong Ha line reported by Reipurth & Aspin (2004) and Briceho et al. 
(2004) also exhibits a P Cygni profile. These profiles are clear indications of mass outflow 
in a wind. The Brackett emission lines, however, lack any absorption component and the 
centroid of the Bra emission is slightly blueshifted, consistent with the findings of Najita 
et al. (1996). This variation in H line profiles with recombination level is often seen in the 
spectra of hot massive stars with strong ionized and accelerating stellar winds, such as the 
B-type supergiants HDE 316285 (Hillier et al. 1998) and P Cygni itself (Najarro, Hillier, & 
Stahl 1997; Lenorzer et al. 2002). HDE 316285 also exhibits permitted lines of Ca II and Na 
I in emission in its NIR spectrum. The estimated outflow velocity from MNO also matches 
the velocities seen in these systems. The similarities between the NIR spectra of these B 
supergiants and that of MNO provide support for the ionized wind model suggested by Nisini 
et al. (2004). Additional support is provided by the B-type optical spectrum of region 'C in 
the reflection nebula presented by Briceno et al. (2004). Presumably, the optical spectrum 
arises from the atmosphere of the disk itself, which may be driving the ionized wind. 

Attempts to use the Pa and Br emission lines to corroborate the extinction estimate 
derived from the ice feature reveal a remarkable discrepancy with the standard Case B 
conditions for H I recombination (Hummer & Storey 1987). The estimated Ay value depends 
strongly on which line ratios are used, ranging from a minimum of 0.5 for the Br7/Bra ratio 
to a maximum of 30 for the Pa<5/Pa7 ratio. The Br line ratios in particular suggest that 
the optically thin Case B assumption is violated for at least some of the low-lying H levels 
and transitions. A similar result has been reported for HH 100 by Nisini et al. (2004). We 



-7- 



expect that some of the Balmer lines may also be optically thick, although clearly a full 
wind model is needed to determine the H level populations, from which the optical depths 
and line profiles can be predicted. High optical depth in the ionized wind would also explain 
why the 3.6, 6, and 20 cm radio emission has not increased during the outburst (Andrews 
et al. 2004). 

Under the assumption that the Br transitions are optically thick, their strengths provide 
an estimate of the size of the emitting regions. Following Nisini et al. (2004), we estimate 
the size of the Br emission region by assuming the line flux is given by blackbody emission 
and adopting a gas temperature of 10 4 K. We find that the Br a emission region has a radius 
~ 8.5i? . From the correlations between mass loss rate and bolometric luminosity, and 
between mass loss rate and Br a luminosity, presented by Nisini et al. (1995), the inferred 
mass loss rate for MNO is of the order of 4 x lO~ 8 M yr _1 , substantially below the typical 
values found for FUor objects (Hartmann & Kenyon 1996). From these values and the 
wind velocity determined from our spectrum, we estimate a density in the emission region 
of the order of 10 9-10 cm~ 3 . For such a dense wind, collisions will play a significant role in 
determining the populations of the low-lying H levels and excited-state photoionization will 
be important in maintaining the ionization equilibrium (see e.g., Simon et al. 1983). A similar 
explanation has been suggested for the H emission lines in the spectrum of the embedded 
B1.5Ve star MWC 297, which may also possess a disk that drives an outflow (Drew et al. 
1997). At these high densities, the forbidden lines normally seen in photoionized regions 
may be collisionally de-excited, which would explain why they are not seen in the optical 
spectrum of MNO (Reipurth & Aspin 2004; Briceno et al. 2004). 

Although the connection between MNO and other eruptive objects (FUors and/or EX- 
ors) is unclear at the moment, the conclusion that MNO is a Class I object seems secure. It 
is possible that the picture we have outlined for MNO may be applicable to other Class I 
sources in L1630, and MNO may therefore represent a natural evolutionary transition from 
the deeply embedded Class objects in nearby areas of the same cloud to the visible Class 
II T Tauri stars that are also present in the vicinity. 

WDV thanks Goran Sandell for many useful discussions and Sean Andrews and Barry 
Rothberg for sharing their results prior to publication. We thank Dave Griep and Bill Golisch 
for their assistance at the telescope. This work was supported by the National Aeronautics 
and Space Administration under Cooperative Agreement no. NCC 5-538 issued through the 
Office of Space Science, Planetary Astronomy Program. 
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Table 1. Line fluxes, EWs, and luminosities in the NIR spectrum of MNO from 2004 

March 09, 06h UT. 



Line 


Wavelength 


Flux a 


EW a 


Luminosity b 




(/iin) 


(W/m 2 ) 


(A) 


(W) 


Ca II 


0.8498 


1.6 x 10~ 17 


-9.8 


6.7 x 10 24 


Ca II 


0.8542 


2.1 x 10~ 17 


-12.4 


8.2 x 10 24 


Ca II 


0.8662 


1.8 x 10~ 17 


-10.1 


6.1 x 10 24 


Pa 5 c 


1.0049 


2.7 x 10~ 18 


-0.6 


2.5 x 10 23 


He I d 


1.0830 




8.9 




Pa 7 C 


1.0938 


2.1 x 10~ 17 


-2.8 


1.2 x 10 24 


OI(?) 


1.1287 


1.6 x 10~ 17 


-1.8 


7.8 x 10 23 


Pa c 


1.2818 


1.7 x 10~ 16 


-10.9 


5.0 x 10 24 


Mgl 


1.5031 


5.4 x 10~ 17 


-2.2 


7.9 x 10 23 


Pa a c 


1.8751 


2.1 x 10" 15 


-52 


1.5 x 10 25 


Br 5 + Ca I 


1.9446 


1.4 x 10- 16 


-3.4 


9.6 x 10 23 


He I d 


2.0581 




0.7 




Br 7 


2.1655 


2.3 x 10~ 16 


-4.6 


1.4 x 10 24 


Na I 


2.2084 


4.8 x 10~ 17 


-0.9 


2.9 x 10 23 


Br a 


4.0512 


6.6 x 10- 16 


-12 


1.9 x 10 24 


CO 


4.67 




65 





a Estimated errors are 10% in both flux and EW. 

b Dereddened luminosity for Ay = 11 and d = 400 pc. 

C P Cygni profile. Flux and EW refer to emission component 
only. 

d P Cygni profile. Flux and EW refer to absorption component 
only. 
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Fig. 1.— (left)The 0.8-5.2//m spectrum of McNeil's Nebula Object obtained with SpeX. 
Prominent emission and absorption features are labeled, (right) Normalized J- and K-band 
spectra of MNO. 



